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The evolution of on-site number fluctuations of ultracold atoms in optical lattices is experimentally 
investigated by monitoring the suppression of spin-changing collisions across the superfluid-Mott 
insulator transition. For low atom numbers, corresponding to an average filling factor close to 
unity, large on-site number fluctuations are necessary for spin-changing collisions to occur. The 
continuous suppression of spin-changing collisions is thus a direct evidence for the emergence of 
number- squeezed states. In the Mott insulator regime, we find that spin-changing collisions are 
suppressed until a threshold atom number, consistent with the number where a Mott plateau with 
doubly-occupied sites is expected to form. 

PACS numbers: 03.75.-b,03.75.Hh 



One of the most fundamental signatures of the Mott 
insulator (MI) transition undergone by ultracold atomic 
ffases in optical lattices [1, J, 1, i, i, 0, 0, B i, E 



11 



121 . Il3| is a drastic change in atom number statistics. In 
a very shallow lattice, ultracold bosons tend to form a 
Bose-Einstein condensate. In this case, a measurement 
of the probability for finding n atoms at a given lat- 
tice site would reveal a characteristic Poisson distribu- 
tion with large on-site fluctuations. However, for deeper 
lattices, the influence of repulsive interactions, which dis- 
favor such fluctuations, becomes increasingly dominant 
and results in the emergence of number-squeezed states 
with suppressed number fluctuations. Above a critical 
lattice depth, the ultracold gas enters the MI regime, 
where the number fluctuations almost vanish. In ex- 
periments so far, interaction-induced number-squeezed 
states were detected through the observation of increased 
phase fluctuations - the canonically conjugate variable to 
number fluctuations 0, [l, through an increased 

timescale for phase diffusion [J]. 

In this Letter, we directly observe the continuous sup- 
pression of number fluctuations when the ultracold sam- 
ple evolves from the superfluid (SF) regime to deep in 
the MI regime. The idea behind our measurement is il- 
lustrated in Fig. [H After producing an ultracold gas in 
an optical lattice, we suddenly increase the lattice in- 
tensity, suppressing tunneling and freezing the number 
distribution. A probe sensitive only to the presence of 
atom pairs at a given lattice site is finally applied. Close 
to unity filling, a non-zero probe signal is obtained only 
if initially large on-site fluctuations produce a non-zero 
fraction of sites with two atoms. While we observe this 
behavior for a gas initially in the SF regime, the probe 
signal is progressively suppressed when approaching the 
Mott transition, indicating increasingly number-squeezed 
states. 

The specific two-particle probe used in this work are 
spin-changing collisions (see [ij] and references therein), 
which convert at each lattice site pairs of spin / = 1 




• * * 



FIG. 1: Illustration of the number statistics measurement. 
Spin changing collisions turn atom pairs initially in the Zee- 
man substate m = (no arrow) to pairs in m = ±1 states (up 
and down arrows). This process happens for sites with n — 2 
(a) or n = 3 (b) atoms. For one atom per site on average, 
whether this occurs depends drastically on the many-body 
correlations. For a Bose-Einstein condensate (c), large on- 
site fluctuations create a finite number of sites with 2 or 3 
atoms, where ±1 pairs can be created. On the contrary, for a 
MI state (d), only isolated atoms are found and no m = ±1 
pairs are created. 



atoms in the m = Zeeman sublevel to pairs with one 
atom in m = +1 and the other in m = — 1. In principle, 
other schemes, e.g. measuring the interaction energy 11 , 
or monitoring atom losses due to Raman photoassocia- 
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16[ or Feshbach resonances [l7|, could be suit- 
able for this measurement. Spin changing collisions are 
appealing because non destructive (see also [18]), and be- 
cause they can be resonantly controlled using the differ- 
ential shift between Zeeman sublevels induced by an off- 
resonant microwave field ^19, 20, 24]. We show that this 
technique allows to measure selectively doubly-occupied 
sites in the optical lattice. 

Our experimental setup has been described in detail in 
[l3|. We first load a degenerate gas of ^^Rb atoms in the 
\F = 1, m = —1) Zeeman sublevel into a combined mag- 
netic trap plus optical lattice potential at an initial lattice 
depth Vq. The intensities are then rapidly increased from 
to Vf = 40£^r within tup = 1 ms (see Fig. [2]). Here 
Er = h? I^My? is the single photon recoil energy, and 
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A = 842 nm the lattice laser wavelength. Immediately af- 
ter this ramp, the magnetic potential is switched off, and 
the cloud is held for 60 ms in order to let the magnetic 
bias field stabilize to its final value B ^ 1.2 G. The atom 
are then prepared in the m = state using microwave 
transfer pulses, and held for a variable time tosc, dur- 
ing which a collisional spin oscillation takes place. This 
coherent evolution is detected experimentally as a re- 
versible exchange between the populations in the m = 
and m = ±1 Zeeman sublevels, measured by absorption 
imaging after 12 ms of free expansion. 
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FIG. 2: Time sequence of the experiment. 

Spin changing interactions in deep optical lattices have 
been described in details in [l^ . We consider spin / = 1 
atoms and assume that tunneling can be neglected, so 
that the lattice sites are isolated from each other. At a 
single lattice site, the spin-changing collisions are criti- 
cally sensitive on the filling n of the well. For sites with 
filling n = or 1, spin-changing collisions cannot oc- 
cur. The first non-trivial cases corresponds to doubly- 
occupied wells. In this case, only two spin states (one 
with both atoms in m = and the other with a single 
m = ±1 pair) are accessible (see Fig.fTJi). Therefore, the 
atom pair undergoes Rabi-like oscillations at the effective 
Rabi frequency Y^of-hT^f- The energy mismatch ("de- 
tuning") between the two states is M2 = Ae + t/g, where 
Ae = e+i + e_i — 2eo corresponds to the difference in Zee- 
man energies Cm- The spin-dependent interaction energy 
Us depends on atomic and lattice parameters [l3], and 
also determines the coupling strength as 1^2 = 2y/2Us. 
Sites with n = 3 atoms behave in a similar way (Fig.[TK), 
however with an energy difference M3 = Ae — Us and a 
coupling strength fts = 2>/6/7s. 

In principle, site occupancies n > 3, whose spin dy- 
namics involve more than one m = ±1 pair, are also 
possible. However, during the hold time thoid indi- 
cated in Fig. [21 those sites can be emptied by three- 
body recombination (3BR) events at an event rate 7^ = 
73Bn(n — l)(n — 2) [21[,with73B ~ 0.5 s~^ for our param- 
eters. Therefore, sites with n > 4 are efficiently removed 
after the wait time. 

In our experiment, we produce large ensembles of 
atoms in the optical lattice with spatially inhomogeneous 
atom number distribution. The inhomogeneity results 
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FIG. 3: Resonance curve of the spin-oscillation amplitude in 
the far-detuned microwave field, measured at a magnetic field 
B ^ 1.2 G and at a fixed hold time thoid = 15.5 ms. Filled 
(hollow) circles denote the population in the Zeeman substate 
m = (resp. in m = il)- 



from an additional trapping potential Ve^t present on 
top of the optical lattice [28]. In the MI regime, this 
potential leads to the formation of fiat Mott plateaux 
with well-defined atom number per site jH, 0, [23]- Also, 
the local fiuctuations have an inhomogeneous distribu- 
tion. Experimentally, we measure the "spin-oscillation 
amplitude" for the entire atomic cloud, i.e. the global 
population ^osc = (A^+i + A^-i)/A^ of the m = ±1 states 
after an evolution time tosc, normalized to the total atom 
number N. This amplitude is related to the probability 
Vn of finding n atoms per lattice site, averaged over the 
cloud spatial profile. 

Let us suppose that we are able to tune the single- 
particle detuning to Ae = —/7s, such that doubly- 
occupied sites are exactly on resonance. Then, neglecting 
sites with n > 4, the oscillation amplitude is obtained by 
summing the contribution from sites with n = 2 and 
n = 3, 
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From Eq. ([T]), we conclude that atom pairs and triplets 
oscillate essentially out of phase. By choosing Ae = —Us 
and t-j^ = 7r/r^2, all doubly-occupied sites are converted 
to m = ±1 pairs, whereas the conversion efficiency for 
triplets is around 3 %. Recording the amplitude of the 
spin oscillations thus allow to probe the distribution of 
atom pairs alone. This is reminiscent of cavity quantum 
electrodynamics 22, [23|, where Fock states of the cav- 
ity field could be discriminated due to different coupling 
strengths to an atomic transition. In particular, choosing 
Ae = Us would allow to measure the fraction of triply- 
occupied sites remaining after three-body decay. 

To achieve full conversion of doubly-occupied sites, it 
is necessary to tune the spin oscillations for doubly oc- 
cupied sites into resonance, i.e. set Ae = —Us. In a 
magnetic field 5, the quadratic Zeeman shift contributes 
a positive amount to Ae. Hence, if Us > (which is 
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FIG. 4: Amplitude of the spin oscillation vs atom number and different lattice depths: Vo = 4, 8, 11, 13, 20, 40 Er (a-f). The 
thin dashed lines show the prediction of a theoretical model, where the atom number distribution is deduced from a mean field 
approach at T = (see text). The thin solid lines show only the fraction of atom pairs calculated from the same model. The 
thick solid lines are guides to the eye. In (f), the vertical lines indicate where Mott plateaus with 2 and 3 atoms per site are 
expected to form. 



the case for ^^Rb), the interaction energy Us leads to a 
residual detuning in zero magnetic field that prevents to 
reach the resonance. For this reason, we introduce a dif- 
ferent technique using the differential level shift induced 
on the individual Zeeman sublevels by a far off-resonant 
microwave field ("AC-Zeeman shift"). With a suitable 
choice of polarization, detuning and power, the detun- 
ing Ae can be tuned at will in the range of interest, and 
allows to compensate the magnetic field contribution to 
Ae plus the interaction term Us- In this work, the mi- 
crowave field is detuned by several hundred MHz to the 
red of any hyperfine resonance to suppress population 
transfer to / = 2. Indeed, no such transfer is observed 
within our experimental sensitivity. 

In Fig. [3l the fraction of atoms found in m = ±1 
is plotted as a function of the microwave power for a 
fixed tosc = 15.5 ms ~ t^^, corresponding to maximum 
conversion. These data were taken for constant initial 
lattice depth (Vb = Vf = 40 Er) and atom number 
{N ^ 2.6 X 10^). For very low microwave powers, the 
spin dynamics is suppressed by the quadratic Zeeman 
detuning (Ae ^ 27r x 207 Hz), much larger than the 
spin-dependent interaction Us ^ 27r x 10.7 Hz. The AC- 
Zeeman shift can compensate for this detuning and for 
the interaction part, inducing a resonance in the num- 
ber of m = ±1 pairs shown in Fig. [3l The oscillation 
amplitude, close to the expected ^osc ^ 0.5 indicates 
that nearly all atom pairs are converted into ±1 pairs, in 
agreement with further experiments discussed in a com- 
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We now turn to the measurement of number statis- 
tics. We choose tosc ^ 15.5 ms and and a dressing field 
tuned to resonance, as in the previous paragraph. At 
a given lattice depth Vb, we have recorded the oscilla- 
tion amplitude in a broad range of atom numbers, from 
about 10^ to a few 10^ [25i,]. The experiment is then 
repeated for various lattice depths, from the SF regime 
{Vo = 4:Er) to deep in the MI regime (Vb = '20 Er and 
40 Er). As shown in Fig.Hl at low lattice depths, the spin 
oscillations occur for any atom number with an am- 
plitude slowly increasing with N. For small atom num- 
ber, the oscillation amplitude is increasingly suppressed 
with increasing lattice depth, and completely vanishes for 
large lattice depths. This qualitative behavior is consis- 
tent with the behavior expected from the Bose-Hubbard 
model S 0, 0, 13, i, El, [H, m . On approaching the 
Mott transition, the ground state adapts to an increased 
interaction energy by reducing its number fluctuations, 
eventually producing an array of one-atom Fock states at 
each site where spin-changing collisions cannot occur. 

Within the MI regime (Fig. [Hi-f), we observe that 
the suppression of spin oscillations persists up to some 
threshold atom number (6.0(3) x 10^ for the data in 
Fig. Hf). This is consistent with the expected formation 
of Mott plateaus with increasing atom number, as the 
cloud expands in the trapping potential. A Mott plateau 
with n atoms per site forms when the cloud radius reaches 
the size Rn where the potential energy V^xt(^n) matches 
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the on-site interaction energy U{n — 1). For a harmonic 
potential with trapping frequency c^ext, this happens at 
a threshold number ]^ ^ ^2 ELi Above 
N2 ~ 47r/3(mcjgx^(i^/2[/)~^/^, a core with two atoms 
per site starts to grow, thus enabling the spin oscilla- 
tions. For the parameters that correspond to Fig. [If 
(^ext = 27r X 80 Hz and Vq = 40 Er)^ we calculate 
A^th ^ 6.8 X 10^, close to the measured value. For even 
higher atom number (corresponding to TVs ~ 3 x 10^), a 
shell of triply-occupied sites start to form, reducing the 
fraction of atoms in the n = 2 shell. This can be seen in 
Fig. [It, where we indeed observe a decrease of the spin 
amplitude above this number. 

In order to compare our experimental results with the 
prediction of the Bose-Hubbard model [5], we solve this 
model numerically within a mean-field approximation at 
zero temperature [sol, [3l|. Accounting for losses dur- 
ing the hold time thoid (wait time plus oscillation time), 
we obtain the distribution Vn- For each filling n and a 
given tosc, we calculate the conversion efficiency r]±i{n) to 
m = ±1 pairs, and obtain the total spin amplitude from 
nri±i{n)Vn- The results of this calculation, 
indicated by the solid line in Fig. [H lie very close to the 
fraction of pairs V2 predicted by the same model (dashed 
line), in agreement with the arguments leading to Eq. [H 
Deep in the MI regime (Fig.[4t-f), the calculations agree 
well with the measurements. For lower lattice depths, al- 
though the qualitative trend is still reproduced, we find 
discrepancies. Near the Mott transition (Fig. [Hi), the 
mean-field calculations predict an amplitude lower than 
observed, a behavior consistent with the study of num- 
ber correlations beyond mean-field reported in [l2[. Be- 
low the transition point (Figs. [4^-c), the model predicts 
an oscillation amplitude higher than observed. Devia- 
tions from the initial distribution may arise in this low 
lattice depth regime if excitations are generated during 
the preparation phase and result in an increased pop- 
ulations in the low-density regions of the cloud, which 
barely participate to the spin oscillations. Such "finite 
temperature" effects have possibly less inffuence in the 
MI regime, where the many-body system is protected by 
an interaction gap. 

In conclusion, we have shown how spin oscillations can 
be used to probe number squeezing in optical lattices via 
the detection of the fraction of atom pairs. Our obser- 
vations confirm the expected scenario: near-Poissonian 
fiuctuations for shallow lattices, strongly suppressed fiuc- 
t nations for deep lattices, and a smooth interpolation in 
between. Moreover, the observed behavior is consistent 
with the expected formation of Mott plateaus, a signa- 
ture of the incompressibility of this system. Our results 
indicate that number squeezing is robust with respect 
to experimental manipulations, such as transfer to the 
purely optical trap. In this sense, they are promising to 
employ those number squeezed states, e.g. in Heisenberg- 
limited atom interferometry [32]. 
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